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Abstract 

We derive the general upper bounds on the mass of the lightest neutralino, as a 
function of the gluino mass, in different supersymmetry breaking models with minimal 
particle content and the standard model gauge group. This includes models with 
gravity mediated supersymmetry breaking (SUGRA) , as well as models with anomaly 
mediated supersymmetry breaking (AMSB). We include the next-to-leading order 
corrections in our evaluation of these bounds. We then expand the mass matrix in 
powers of Mz/n, and find the upper bound on the mass of the lightest neutralino 
from this expansion. By scanning over all of the parameter space, we find that the 
bound we have obtained can be saturated. We compare the general upper bound on 
the lightest neutralino mass to the upper bound that is obtained when the radiative 
electroweak symmetry breaking scenario is assumed. 



PACS number(s): 12.60.Jv, 14.80.Ly 



1 Introduction 

It is widely expected that at least some supersymmetric particles will be produced at the 
large hadron collider (LHC) that is starting operation in a few years' time. However, 
most of these supersymmetric particles will not be detected as such, since they will decay 
into the particles of the Standard Model (SM), or to the lightest supersymmetric particle 
(LSP), which is stable as long as the R-parity is conserved. Thus, the experimental study 
of supersymmetry involves the study of cascade decays of the supersymmetric particles 
to the LSP and the reconstruction of the subsequent decay chains. The LSP in a large 
class of supersymmetry breaking models is the lightest neutralino, which has thus been a 
subject of intense study for a long time PQ El El III El E] • A stable lightest neutralino is also 



an excellent candidate for dark matter [5]. As such it is important to have information on 
the mass of the lightest neutralino state. 

In view of this, the properties of the lightest neutralino and also heavier neutralinos and 
charginos, which often appear in the cascade decays, are of considerable importance. In 
the minimal version 1 of the supersymmetric extension of the Standard Model at least two 
Higgs doublets Hi and H2 with hypercharge (Y) having values —1 and +1, respectively, are 
required. The fermionic partners of these Higgs doublets mix with the fermionic partners 
of the gauge bosons to produce four neutralino states x®,i = 1,2,3,4, and two chargino 
states xt > * = 1) 2. 

The neutralino mass matrix A4 depends on the ratio of the vacuum expectation val- 
ues (VEVs) of the two Higgs doublets denoted by tan/3 = V2/V1, where v\ = (H®) and 
i>2 = (H®) are the vacuum expectation values of the two Higgs doublets with opposite 
hypercharge, the supersymmetry breaking U(l)y and SU{2)i gaugino masses M\ and 
M2, and the supersymmetry conserving Higgs (ino) mixing parameter /U. The mass matrix 
is symmetric, but not necessarily real. The mass parameters can have arbitrary complex 
phases, as can also the Higgs boson VEVs. However, all of these are not actually indepen- 
dent - one can choose the two nontrivial phases to be in Mi and fj,. The electric dipole 
moments strictly constrain the phases in supersymmetric (SUSY) models. However, these 
bounds are for products of the phases. Thus, if there are cancellations between phases, a 
single phase can be larger than the limits for the product [7j. 

Recently, it has been demonstrated that at the linear collider one can determine the 
above parameters of the neutralino and chargino sectors from the masses of charginos 
and three lightest neutralinos, or alternatively from two lightest neutralinos and the cross 
section e + e~ — > X1X2 [DEI- The linear collider is likely to be available several years after 
the completion of the LHC, and thus all the information that is available now or can be 
obtained at the LHC will be very valuable. 

In this paper we obtain the theoretical upper bound on the mass of the lightest neu- 
tralino state in the most commonly studied supersymmetry breaking models. These in- 
clude the gravity mediated supersymmetry breaking model and the anomaly mediated 
supersymmetry breaking model, with minimal particle content. 2 In a general model with 
an arbitrary particle content, an upper bound for the lightest neutralino mass was calcu- 
lated in For specific supersymmetry breaking scenarios we can give a more accurate 
bound. In Section 2 we obtain the general upper bound on the mass of the lightest neu- 
tralino in the minimal version of the supersymmetric standard model. We then evaluate 
this upper bound for the two most popular supersymmetry breaking models, namely, the 

x By minimal version we here mean the model with the minimal particle content and the Standard 
Model gauge group 

2 In the gauge mediated supersymmetry breaking (GMSB) models the lightest neutralino is not the 
lightest supersymmetric particle. However, in many models it is the next-to-lightest particle. Here we will 
also comment on the upper bound in the GMSB models. 



gravity mediated super symmetry breaking (SUGRA) models and the anomaly mediated 
supersymmetry breaking (AMSB) models. We include the next-to-leading order correc- 
tions in the numerical evaluation of this upper bound. 

In Section 3, we then study the expansion of the neutralino mass matrix in powers 
of Mz/fJ>, which at second order is accurate to 1% for large values of /i. The extremum 
(maximum) of the lightest neutralino mass gives an upper bound from this expansion. 
This value is lower than the upper bound obtained directly from the mass matrix. By 
a numerical scan over real and complex parameter values we confirm that this bound is 
accurately saturated in the supersymmetric models that we study in this paper. We com- 
pare the general upper bound with the largest value of the lightest neutralino mass that is 
obtained when radiative electroweak symmetry breaking is assumed (and the \x parameter 
determined from the radiative electroweak symmetry breaking). For this purpose we have 
used the numerical program SOFTSUSY 9 . In Section 4 we present our conclusions. 



2 The general upper bound on the mass of the lightest neu- 
tralino 

We start by recalling the neutralino mass matrix in supersymmetric models in the basis 
jjP. = (-i\>, -a 3 , rP X Hl , iP 2 h 2 ), 3 = 1, 2, 3, 4, (1) 
which can be written as ^U] 



M 



f Mi -M z cos (3 sin 9 W M z sin (3 sin 9 W \ 

M2 Mz cos (3 cos 9\y —Mz sin f3 cos 9w 

— Mz cos /3sin 9w Mz cos (3 cos 9w —\i 

1 Mz sin (3 sin 9w —Mz sin (3 cos 9w — /i 1 



, (2) 



where A' and A3 are the two-component gaugino states corresponding to the U(1)y and 
the third component of the SU(2)l gauge groups, respectively, and V^Zi > are the two- 
component Higgsino states. Furthermore, g 1 and g are the gauge couplings associated 
with the U{1)y and the SU(2)l gauge groups, respectively, with tan#vy = g'/g, and 
M| = (g 2 + g' 2 ){v\ + «2)/2. Assuming CP conservation, this mass matrix is real. We shall 
denote the eigenstates of the neutralino mass matrix by XvX2^X3iXi labeled in order of 
increasing mass. Since some of the neutralino masses resulting from diagonalization of 
the mass matrix can be negative, we shall consider the squared mass matrix Ad^Ai. An 
upper bound on the squared mass of the lightest neutralino Xi can be obtained by using 
the fact that the smallest eigenvalue of Ai^Ai is smaller than the smallest eigenvalue of 
its upper left 2x2 submatrix 



/ M\ + Mf sin 2 9 W -Mf sin 9 W cos 9 W N 
v -M| sin 9 W cos W M| + Mf cos 2 6>h/ y 



(3) 



thereby resulting in the upper bound 

M 2 o < ~ f Mf + Ml + Mf - ^ (Mf - M|) 2 + M| - 2(Mf - M 2 2 )Mf cos 20^ . (4) 

We emphasize that the upper bound @ is independent of the super symmetry conserving 
parameter and also independent of tan [3, but depends on the supersymmetry breaking 
gaugino mass parameters M\ and M2. Despite this dependence on the unknown super- 
symmetry breaking parameters, we will show that Eq. Q leads to a useful bound on 

M x\- 

2.1 Gravity mediated supersymmetry breaking 

In the gravity mediated minimal supersymmetric standard model, the soft gaugino masses 
Mj satisfy the renormalization group equations (RGEs) (IM3I = m g , the gluino mass) 

l ^ d J^ = 2hl M t gl h= (f,l,-3) (5) 

at the leading order. Here g\ = |</, <?2 = g, and 53 is the SU(3)c gauge coupling. The 
RGEs © imply that the soft supersymmetry breaking gaugino masses scale like gauge 
couplings: 

M X {M Z ) = M 2 (M Z ) = M 3 (M Z ) 

«i(M z ) a 2 (M z ) a 3 {M z ) ' [ ' 

where «j = gf /4tt, i = 1,2, 3. 

The relation © reduces the three gaugino mass parameters to one, which we take to 
be the gluino mass rag. The other gaugino mass parameters are then determined through 



5 a 

3a3 cos 2 9w 
a 

Q?3 sin^ 9w 

where we have used the value of various couplings at the Z° mass 



M i{M z ) = — m g ~ 0.14 ms, (7) 

M 2(M Z ) = ™_2 * m 9 ~ °- 28 m 9' ( 8 ) 



a~ l (M z ) = 127.9, sin 2 9 W = 0.23, a 3 (M z ) = 0.12. (9) 

Using Eqs. (0) and © in Eq. (J3J), we get the upper bound on the mass of the lightest 

neutralino. 3 For a gluino mass of 200 GeV, the upper bound @ for the lightest neutralino 

3 In the GMSB models, one gets the same relations JSJ at the messenger scale, since M a oc a a . Thus, 
the upper bound obtained in the SUGRA model can be applied in the GMSB model as well. 
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Figure 1: The upper bound on the mass of the lightest neutralino in the SUGRA and 
AMSB models. The tree level results are given by dashed lines and the next-to-leading 
order results by the solid lines. 



mass is about 35 GeV. Similarly, for a gluino mass of 1 TeV, the upper bound Q becomes 
186 GeV. 

We have plotted the upper bound (jlj) on the mass of the lightest neutralino in Fig. ^ 
The almost straight dashed line corresponds to the SUGRA model at the tree level. From 
Fig. n we observe that m x o < 186 GeV for m g < 1 TeV. 

We now include next-to-leading order (NLO) corrections coming from 03 and from the 
top-quark Yukawa coupling at (= h^/Air) two- loop contributions to the beta functions and 
logarithmically enhanced weak threshold corrections. In this approximation, one finds jllj 

f a T O 2 m? rr 2 

Mi(Q) \ 1 + 



M NLO 



in cos 2 6w 



Q 2 ™J m l 

-211n^ + 11 In— I + 91n— | 

M\ Ml Ml 



+ ln ^2 +TT si * 2 P 



2fi 



m. 



In 



2 



+ 



2«3 13«t 



3ir 66ir 



(10) 
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M 2 (Q) I 1 + 
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8tt sin 2 1 
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M| Ml Ml 



+ l n "tto + tt" sin 2/3 — ^ — 

M 2 Mi H [i 2 -m\ 
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(11) 



F(x) 




M 3 (Q) l + £ 



1 + 2x + 2x(2 - x) lnx + 2(1 - x)^ ln|l 



(12) 
(13) 



Here Mi{Q) are the leading order results given by Eq. @. Notice that the next-to-leading 
order corrections are of the same form in all models. It is only the leading order Mi(Q) 
that are different for different models. 

In order to calculate the next-to-leading order upper bound on the lightest neutralino 
mass, we need to know Al^ LO and M^ 1 " . As with the leading order result, we express 
M^ LO and M^ LO as a function of M^ LO (the NLO physical gluino mass), using Eqs. (|10|) 
- lfT5|) and substitute it in Eq. (j3J). We plot the NLO corrected upper bound on the lightest 
neutralino mass as a function of the gluino mass as a solid curve in Fig. ^ As input we 
have used here tan/3 = 10, m t {pole) = 17 '4 GeV, m = 300 GeV, A = 1 TeV, p = -460 
GeV, and Q = 890 GeV. Since dependence on these model parameters appears only at 
the loop level, the upper bound is not very sensitive to these parameters. We note that 
the NLO corrections increase the upper bound from its tree level result by only a few 
GeV for a wide range of the gluino mass. Indeed, we find that at the NLO the upper 
bound on the mass of the lightest neutralino is m^o < 194 GeV for rrig < 1 TeV. In ^2] 
full one- loop corrections to sparticle masses were calculated. The loop corrections to the 
lightest neutralino mass can be typically 10% for m v o < 40 GeV. However, for m v o = 100 
GeV the full one loop corrections to the lightest neutralino mass are < 5%. We note that 
the experimental lower bound ^3] on the lightest neutralino mass, valid for any tan (5 and 
mo, is m x o > 37 GeV. 

2.2 The anomaly mediated supersymmetry breaking 

The anomaly induced soft terms are always present in a broken supergravity theory, re- 
gardless of the specific form of the couplings between the hidden and observable sectors. 
They are linked to the existence of the superconformal anomaly. Indeed, they explicitly 
arise when one tries to eliminate from the relevant Lagrangian the supersymmetry break- 
ing auxiliary background field by making a suitable Weyl rescaling of the superfields in 
the observable sector. 

The soft terms in the anomaly mediated supersymmetry breaking models are especially 
interesting because they are invariant under the renormalization group transformations. 
The phenomenological appeal of the soft terms in AMSB resides precisely in this crucial 
property. In particular, it implies a large degree of predictivity, since all the soft terms 
can be computed from the known low-energy SM parameters and a single mass scale 
777.3/2- Also, it leads to robust predictions, since the RG invariance guarantees complete 
insensitivity of the soft terms to the ultraviolet physics. As demonstrated with specific 
examples in Ref. |14| . heavy states do not affect the low-energy parameters, since their 
effects in the beta functions and threshold corrections exactly compensate each other. 
This means that the gaugino mass prediction 

M A = ^m 3/2 (14) 



is valid irrespective of the grand unified theory gauge group in which the SM may or 
may not be embedded. A unique feature of the anomaly mediated supersymmetry is the 
gaugino mass hierarchy implied by Eq. (|14|). At the leading order, we thus have 



M i(Q) = A Uai ?l m 3/2 ^8.9x HT 3 m 3/2 , (15) 
4-7T cos z Ow 



M 2 (Q) = a{Q 2 \ m 3/2 ~ 2.7 x Hr 3 m 3/2 , (16) 
4-7T sin 0\y 



M 3 (Q) = _^M m3/2 ~ _ 2 .8 x 10- 2 m 3/2 , (17) 

at the scale Mz- Using Eqs. ()15jl - (|17|1 in Eq. @, we obtain the leading order result 
for the upper bound on the lightest neutralino mass in the minimal AMSB model. We 
have plotted this upper limit as the upper dashed curve in the Fig. ^ It is interesting to 
note that there is a kink in this dashed curve around rrig ~ 210 GeV. This is due to the 
competition between the diagonal terms in the 2x2 submatrix 0. The term containing 
M\ is smaller, when the gluino mass is small, but with the increasing gluino mass the 
term with M 2 becomes smaller around 210 GeV. This is because the Wino triplet mass 
parameter is always smaller than the Bino mass parameter in the AMSB type model, 
contrary to the SUGRA or GMSB type models where the Bino mass parameter is smaller 
than Mi. 

In the next-to-leading order corrections to the lightest neutralino mass in AMSB mod- 
els, the complete sparticle spectrum becomes important. Unfortunately, it turns out that 
the pure scalar mass squared anomaly contribution for the sleptons is negative In 
order to avoid this problem we need to consider other positive soft contributions to the 
spectrum. This can arise in a number of ways, but most of the solutions will spoil the 
RG invariance of the soft terms and the consequent ultraviolet insensitivity. Nevertheless, 
there are various options to cure this problem without reintroducing the flavor problem 

nsmnnmnHi. 

The necessary cure for the slepton masses may also completely upset the mass relations 
for the other particles (as in the case of the model of Ref. 16 J. However, here we will 
simply parameterize the new positive contributions to the squared sfermion masses with 
a common mass parameter Wq, assuming that the extra terms do not reintroduce the 
supersymmetric flavor problem. The low-energy soft supersymmetry breaking parameters 
for the scalars and the trilinear couplings are then obtained from 

"4 

Ay = -^m 3/2 , (19) 

respectively. Using Eqs. (|10|) - ()13f) . we obtain for the anomaly mediated supersymmetry 
breaking models the next-to-leading order results for the gaugino mass parameters as 

m nlo = i.06 Mi (Q), (20) 



m Q = iy s + ^J ra3 / 2 + m °' (18) 



M^ LO = 1.28 M 2 (Q), (21) 
M^ LO = 0.9M 3 (Q), (22) 

where the Mi(Q),i = 1,2,3 (the leading order result), is given in Eqs. I)15 [) -(|17j ) . Here 
we have used as input tan/3 = 10, m t (pole) = 174 GeV, m 3 / 2 = 35 TeV, tuq = 600 
GeV, \i = —600 GeV, and Q = 958 GeV. The Higgsino corrections to M\ and M 2 are 
proportional to [x/Mip and can become very important in models with large fi, as discussed 
in Ref. [H]. 

In Fig. n we have plotted the next-to-leading order upper bound on the mass of the 
lightest neutralino in anomaly mediated supersymmetry breaking models. The NLO re- 
sult, obtained using Eqs. (|20|) - (J22|). is shown as a solid line. The NLO corrections are 
significant, of the order of 20%. The larger NLO correction in the AMSB model as com- 
pared to the SUGRA model is due to the fact that the 03 corrections for the M 2 mass 
parameter are larger than for the M\ parameter. For mg < 1 TeV, the upper bound on 
the lightest gluino mass is 167 GeV, which is considerably less than in the SUGRA case. 

3 Lightest neutralino mass bound from the structure of the 
mass matrix 

We can also obtain information on the neutralino masses by studying the expansion of the 
neutralino mass matrix in terms of the parameter Mz/n- This expansion can be obtained 
most conveniently by using the basis (7, Z°, H®, H®). In this basis the mass matrix is 
given by 



/ Mm 2 + Mo.sL fMo-MiW.sw \ 



M 



M X (? W + M 2 s 2 w 


(M 2 - Mx)c w s w 








(Af a - Mi)c w s w 


M lS 2 w + M 2 c 2 w 


M z 








M z 












-A*c 2 /3 


/J-S2/3 



(23) 



\ u U -/UC 2/3 ns 2 p J 

Here we have used the abbreviations S2fs = sin 2(3, 02/3 = cos 2/3, s\y = sin 2 6w and 
c\y = cos 2 9yy- Let us start by supposing, as before, that all the mass parameters are real. 
The mass matrix is then real and symmetric. 4 The neutralino mass matrix A4 can be cast 
into a form whereby the gaugino and Higgsino mass parameters are only at the diagonal 
positions by a similarity transformation with a matrix A, 

M = A T MA, (24) 



4 We note that in the specific models that we have been considering, SUGRA and AMSB, the phases of 
Mi and M2 are the same (see Eqs. JSJ and 1141 '), So if M2 is real, then Mi is also real. On the other hand, 
the n parameter is in general complex. Complex parameters would imply a non-Hermitian mass matrix, 
giving generally complex eigenvalues. Such a situation can be handled by considering the eigenvalues of 
the matrix M^M. 
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Figure 2: Mass of the lightest neutralino Xi as a function of /i. Solid lines correspond to 
the numerical results from the mass matrix and dashed lines to the second order expansion 
from Eq. q2b|) . The upper two lines represent masses in SUGRA models and the lower two 
in AMSB models. Here tan /3 = 10 and m- q = 900 GeV. 



where 



A 



\ 



(25) 



cw $w 

-sy/ cw 

cos(vr/4-/?) sin(7r/4-/3) 

y -sin(vr/4-/5) cos(vr/4-/5) y 

The mass matrix can then be diagonalized by using perturbation theory. In the SUGRA 
model, for the mass of the lightest neutralino we get, up to terms of 0(Mz / fi) 2 , 



Mi 



M 2 s 2 

m z s w 



sin 2P - ( MhlyMx + *f? 8 w<%r gin 2 2fj 



li \ ~ ■■ M 2 -M 1 

Similarly, for the second lightest neutralino X2 we obtain 



1 

7 2 



(26) 



A2 



Mo 



M 2 r 2 
1V1 Z C W 



sin 2(3 



(27) 



If instead we were considering the AMSB model, Eq. (|27|) would represent the mass of 
the lightest neutralino x?) an d Eq. H26|) would give the formula for the mass of the second 
lightest neutralino. The dependence of the lightest neutralino mass on the specific SUSY 
breaking scenario is due to the fact that the ordering of the gaugino mass parameters is 
model dependent (for AMSB models M 2 < Mi, whereas for SUGRA models Mi < M 2 ). 

In Fig. |21 we plot the mass of the lightest neutralino obtained from the expansion of 
the mass matrix in (M^//x) together with the exact results obtained from the numerical 
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Figure 3: The upper limit on m v o as a function of m a . Solid lines represent masses in 
SUGRA model and short dashed lines in AMSB model. The lower curve in each case 
corresponds to the upper limit obtained from the expansion in (M///t), and the upper 
curve corresponds to the general upper limit obtained from the mass matrix. The long 
dashed curves correspond to the case when the radiative electroweak symmetry breaking 
scenario is implemented. 




evaluation of the lightest neutralino mass from the mass matrix. The results for the other 
neutralinos are very similar in accuracy. The second order tree level expansion is generally 
better than 1% for \u.\ > 450 GeV (nig < 1600 GeV and tan/3 = 10), with the exception 
of small gluino mass, when m v o is very small, thus giving a larger relative error. For our 
purpose it is sufficient to calculate the expansion up to second order in (Mz / £t) • 

Due to the simple functional form of Eqs. Q26JI and (|27|) the extremal values of the 
masses with respect to \i are easily calculated. These functions have only one extremum 
(maximum), which is given (within the limits of validity of the expansions) for the values 
of n 

( M Mf sin 2/3 s 2 w c 2 w \ 
"-- 2 {^ + M-M t )> (28) 

and the maximum mass is then given by the upper bound 

„ r 1 Mlt 2 sin 2 2(3 . . 

m v o<M + - f-, , , M 2na , 29) 

llV1 "I- M—M 

where M = mm(M 1 ,M 2 ),M = max(M 1 ,M 2 ),t = if M = Mi (SUGRA), and t = cfy 
if M = M 2 (AMSB). In Fig. |31 we plot the upper limit on the lightest neutralino mass 
obtained from Eq. (|29jl as a function of nig, for both the SUGRA and AMSB models. We 
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Figure 4: The upper limits of m x o and m x o as a function of nig. Solid lines represent 
masses in SUGRA model and dashed in AMSB model. 

also plot the upper limit obtained from Eq. (@J) in the same figure. The SUGRA results 
are represented as solid lines and the AMSB results as short dashed lines. The lower curve 
of each set corresponds to the upper bound obtained from the expansion in (Mz/n) 2 , and 
the upper curve corresponds to the upper bound obtained from Eq. (jlj). These results for 
SUGRA and AMSB in Fig. |3] are NLO results. We have plotted the results for the value 
of tan (3 = 10. In order to verify the accuracy of these results, we made an extensive scan 
over the parameter space, using both real and complex values of the /U parameter. The 
highest mass obtained from this corresponds extremely well to the upper limit obtained 
from the expansion in (Mz / n) ■ 

We have also made a scan over the parameter space using the SOFTSUSY program 
0, in which the phenomenon of radiative electroweak symmetry breaking (REWSB) is 
implemented. Thus, the fi value in this program is given by the REWSB condition. The 
resulting spectrum includes one- and dominant two-loop corrections. The maximum mass 
obtained for the lightest neutralino is plotted in Fig.[3]as a function of rrig with long-dashed 
lines. The upper long-dashed line corresponds to the SUGRA model and the lower one to 
the AMSB model. One can see that with radiative electroweak symmetry breaking, the 
m x o in the AMSB model is close to the maximum mass obtained from the expansion in 
(Mz/n), while in the SUGRA model with REWSB the m^o obtained is clearly lower than 
the maximum value from the expansion, indicating that [J-rewsb for the SUGRA model 
is not close to the value obtained from Eq. (|28j). 

As in the case for Xi we can search for the upper bound on the mass of the second 
lightest neutralino X2- For light gluinos (lighter than O(60) GeV) the extremum in the 
mass for \2 ls a minimum due to a sign change in the expansion, but for experimentally 




allowed masses the extremum for m^o is a maximum. In Fig. 0] we have plotted the upper 
limits for both the lightest and second lightest neutralino obtained from the expansion 
in {Mzl\x). The solid lines correspond to the SUGRA model, while the dashed lines 
correspond to the AMSB model. For m~ g < 1 TeV, the NLO upper bounds for the second 
lightest neutralino are 440 GeV for the SUGRA case and 419 GeV for the AMSB case. 

4 Conclusions 

In this paper we have studied the neutralino mass matrix for the minimal supersymmet- 
ric model with the aim of obtaining an upper bound on the mass of lightest neutralino. 
Knowledge of the mass of the lightest neutralino is of crucial importance for the super- 
symmetric phenomenology. We have shown that a general limit, valid for arbitrary values 
of parameters, can be obtained from the mass matrix. Even though such a bound depends 
on the supersymmetry breaking parameters M\ and M2, it nevertheless leads to a signifi- 
cant numerical bound on the lightest neutralino mass in the SUGRA and AMSB models. 
We have also obtained an upper bound on the lightest neutralino mass by expanding the 
neutralino mass matrix in terms of the parameter Mz/fJ-- We see that the upper limit from 
this expansion is considerably lower for the AMSB model than for the SUGRA model for 
similar m g . From this analysis we conclude that the upper bound on the mass of lightest 
neutralino is m v o < 200 GeV for mg < 1 TeV. 

In Fig.EJwe have three separate regions for the upper bound on the mass of the lightest 
neutralino: one which is valid in both SUGRA and AMSB case, one which is valid only 
in one of the models, and a third one which is not available for any of the models that we 
have studied. 
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